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The in-core self-powered neutron detector (SPND) acts as a key measuring device for the monitoring of 
parameters and evaluation of the operating conditions of nuclear reactors. Prompt detection and tolerance of 
faulty SPNDs are indispensable for reliable reactor management. To completely extract the correlated state 
information of SPNDs, we constructed a twin model based on a generalized regression neural network (GRNN) 
that represents the common relationships among overall signals. Faulty SPNDs were determined because of the 
functional concordance of the twin model and real monitoring systems, which calculated the error probability 
distribution between the model outputs and real values. Fault detection follows a tolerance phase to reinforce the 
stability of the twin model in the case of massive failures. A weighted K-nearest neighbor model was employed 
to reasonably reconstruct the values of the faulty signals and guarantee data purity. The experimental evaluation 
of the proposed method showed promising results, with excellent output consistency and high detection accuracy 
for both single- and multiple-point faulty SPNDs. For unexpected excessive failures, the proposed tolerance 
approach can efficiently repair fault behaviors and enhance the prediction performance of the twin model. 
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I. INTRODUCTION 


Neutron flux is an important variable in nuclear reactors. 
Monitoring the variation and distribution of the neutron flux 
is essential to maintain the power stability of nuclear reac- 
tors [1]. At present, self-powered neutron detectors (SP- 
NDs) are widely employed in core neutron flux measurement 
systems (CNFMS) of nuclear power plants (NPPs) to accu- 
rately measure neutron flux and provide highly reliable three- 
dimensional (3D) power distribution information. However, 
the risk of failure has increased with the increasing scale and 
complexity of the control of modern nuclear reactors[2]. 


Faulty SPNDs that either completely or partially fail pro- 
vide incorrect monitoring information, which may negatively 
affect both the simple and more advanced functionalities of 
the system, resulting in degradation of the overall system 
performance and an increased risk level [3]. Techniques to 
address these problems can be classified as hardware redun- 
dancy methods, model-based methods, and data-driven meth- 
ods [4, 5]. Generally, hardware redundancy measures in 
which more than three SPNDs are installed to observe the 
neutron flux within a range of space are employed for NPPs 
to improve the reliability of CNFMS. Assuming that any one 
SPND in a neutron measurement channel fails, the additional 
SPNDs still function and maintain high-accuracy measure- 
ments. Although this measure can prevent occasional inter- 
mittent failures from negatively affecting the system, the in- 
ability to implement fault tolerance can lead to false percep- 
tions of the performance of the CNFMS when redundant SP- 
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NDs fail simultaneously. Moreover, redundant SPNDs con- 
tinue to incur high installation and maintenance costs, which 
adversely affect the economics of NPPs [6]. 


Popular model-based methods (such as Kalman filter- 
ing [7], mixed Kullback—Leibler divergence and exponen- 
tial weighted moving average [8], high-gain observer [9], the 
Monte Carlo method [10], and extended state observer [11]) 
handle fault detection using simple mathematical models that 
formulate fault signatures. The results of these methods for 
specific nonlinear systems have been developed under various 
restrictive assumptions. Recently, data-driven methods [12— 
14] have been used to perform fault detection for ker equip- 
ment in NPPs; the main principle is to establish an object 
model based on data analyses and realize accurate detection 
using model outputs with constraints of the evaluation cri- 
teria [15]. Peng et al. [3] constructed mathematical mod- 
els for various detectors using principal component analysis 
(PCA) and achieved fault detection and isolation for SPNDs 
through the square prediction error of linear models and the 
detector validity index based on the reconstruction. Yellapu 
et al. [16] developed a method based on multiscale PCA with 
wavelet transform to reduce the modeling cost and improve 
the sensor fault diagnosis results by calculating the wavelet 
approximation coefficient. An online multiscale data rec- 
onciliation scheme for detecting and isolating sensor faults 
in advanced heavy-water reactors was proposed in [17]; the 
scheme achieved high accuracy under different scenarios. Li 
et al. [18] employed two different fault identification meth- 
ods to locate faults more accurately. One was an improved 
weighted contribution analysis method based on the tradi- 
tional contribution analysis of sensors to Q statistics. The 
second method was based on the sensor validity index ob- 
tained using the iterative reconstruction method. However, 
static linear models are typically unable to detect long-term 
SPND faults. To improve the detection efficiency of practical 


time-varying faults, Chatterjee et al. [19] proposed the use 
of instantaneous cluster statistics to normalize the measure- 
ments of each SPND in clusters and update the PCA model 
using normalized values. This approach resulted in lower 
false alarm rates and higher detection rates for real-time fault 
changes than traditional static models. Experimental results 
show that these approaches are useful in identifying faulty 
instruments. Yu et al. [20] proposed improvements to the 
traditional PCA model through a new corrected reconstruc- 
tion algorithm to reconstruct the principal component and the 
residual space. A cyclic PCA monitoring model was estab- 
lished to accurately detect different types of faults and recon- 
struct fault data. Nageswara et al. [21] performed information 
fusion by combining the ensemble of trees and the support 
vector machine (SVM) algorithm to evaluate the calculation 
error of multiple sensors and the influence of complementary 
and redundant sensors. Most studies have failed to focus on 
the inherent correlation among all SPND signals in the over- 
all detector assemblies. For the dozens of SPNDs existing 
in neutron flux measurement channels, challenges in infor- 
mation integration and correlation analysis remain because of 
difficulties in handling large amounts of data simultaneously. 


As a new enabling technology, a digital twin (DT) can act 
as a mirror of the real world by providing an integrated en- 
vironment for simulating, decision-making, and optimizing 
physical system operations [22, 23]. Because of the power- 
ful computing capabilities and cognitive intelligence of DT, 
developing more refined and scalable models for fault de- 
tection has become possible. Lin et al. [24] developed a 
nearly autonomous management and control (NAMAC) sys- 
tem for advanced nuclear reactors with DT technology. Cai 
et al. [25] proposed analytical techniques based on data and 
information fusion for modeling and developing DT virtual 
machine tools. DT applications in the context of NPPs have 
proliferated recently [26-28]. The development of DT com- 
bined with advanced technologies for detection, control, and 
optimization can significantly improve system performance, 
reliability, availability, maintainability, and operational flexi- 
bility. Therefore, efficient in-core SPND fault detection and 
the maintenance of monitoring systems using twin technolo- 
gies are our major objectives. To continue such interesting 
exploration of research on SPND signals, a twin model us- 
ing DT technology was constructed for the parameter analy- 
sis of an in-core nuclear reactor with large monitoring data 
and complex and changeable operating conditions. In this 
case, the established twin model can extract rich values from 
isolated time-varying data without disturbing the equipment 
on the real physical layer and can simulate the real-time state 
and dynamic characteristics of SPND entities through interac- 
tive data, overcoming the problem that traditional mathemat- 
ical models [29] cannot effectively deduce the fault state of 
multidimensional signals in real time. In addition, several re- 
searchers have devoted their attention to state analysis and to 
the design of DT models of key nuclear reactor components. 
Hu et al. [30] comprehensively described the research status 
and development directions of DT technology in the field of 
advanced nuclear energy, proposed a multidimensional eval- 


uation digital model suitable for application in nuclear reac- 
tors, and preliminarily established the fault diagnosis process. 
Wang et al. [31] proposed a DT system of the out-of-core 
detector assisted with an installation robot to perform real- 
time visualization monitoring of the detector installation and 
replacement process. Gong et al. [32] combined reduced- 
order models with machine learning to create physics-based 
DTs using parameters input in real time to rapidly reconstruct 
the neutron field in the core. Cancemi et al. [33] generated 
primary nuclear components through numerical simulation of 
different plant conditions, which may support the predictive 
maintenance optimization based on plant condition and the 
development of a DT model for improving plant safety and 
availability. 


Because the implementation of the twin model signifi- 
cantly depends on the information transferred from valid data 
[34], the frequent and unexpected fault behaviors of SPNDs 
only provide unrepresentative data, leading to serious distor- 
tions in the outcomes and even the unprecedented collapse of 
the twin model. The decision made after calculations based 
on a large amount of incorrect information fails to realize the 
entire purpose of maintaining normal system function. Most 
recent research methods have focused on fault detection, and 
only a limited amount of research has been conducted to si- 
multaneously achieve fault tolerance in critical detectors at 
NPPs. This is one of the principal motivations for this work. 
Designing a reasonable fault tolerance strategy that would 
work uninterruptedly with sufficient information-processing 
capabilities is practically challenging. Li et al. [35] sug- 
gested an active fault tolerance control method based on the 
improved BP neural network, which controls fault sensors by 
reconstruction. Kim et al. [36] developed a method of eval- 
uating the fault detection coverage of a fault-tolerant tech- 
nique using a fault injection experiment in a safety-critical 
digital I&C system for NPPs. They also proposed a proba- 
bilistic safety assessment model to observe the effect of the 
fault detection coverage of fault-tolerant techniques. To en- 
hance the stability of a modular high-temperature gas-cooled 
reactor system and the control rod drive mechanism, Hui et 
al. [37] proposed an adaptive fault-tolerant control scheme 
based on a radial basis function neural network, which has a 
higher load tracking accuracy and better fault tolerance to sys- 
tems. Li et al. [38] proposed an active fault-tolerant control 
scheme based on the deep Q-network algorithm of reinforce- 
ment learning to maintain the stability of the once-through 
steam generator control system. Rangegowda et al. [39] in- 
troduced a fault tolerance control framework to arrest the per- 
formance degradation of conventional controllers in the pres- 
ence of sensor bias. Clearly, appropriate fault tolerance op- 
erations are required to handle the incorrect information pro- 
vided by failed equipment, which can ensure data purity and 
no interference for the normal performance of the twin model. 
Therefore, introducing an appropriate fault-tolerance strategy 
to SPND fault detection is advantageous to maintain a healthy 
interaction of data information and pave the way for reliable 
outputs of the twin model in the long term. 


In this study, we propose an effective fault detection and 


tolerance approach for in-core SPNDs based on a twin model. 
The SPNDs are uniformly distributed in the reactor core; a 
generalized regression neural network (GRNN) [40] is em- 
ployed to construct the twin model, which is consistent with 
the real system and represents the common relationship be- 
tween the overall signals. In this manner, the state correlation 
between the SPND signals is completely considered, and the 
output characteristics of specific SPND individuals can be de- 
scribed using the joint feature information of the surrounding 
SPNDs. Then, state analysis and fault detection of the SPNDs 
were realized by analyzing the probability distribution of er- 
rors between the outputs of the twin model and the real value. 
To achieve fault tolerance for unexpected faulty SPNDs, we 
useda weighted k-nearest neighbor (WKNN) [41] to recover 
and reconstruct faulty signals. Through troubleshooting, data 
substitution, and data verification of faulty SPNDs, the ac- 
curacy and rationality of the detection results can be signifi- 
cantly improved. Compared with the traditional single-signal 
analysis model, the proposed approach can detect multi-point 
faults more simply and efficiently. 


The remainder of this paper is organized as follows: Sec- 
tion 2 provides a brief description of the measurement chan- 
nel distribution of the neutron flux and the composition of 
the SPND. Section 3 describes the construction scheme of 
the twin model and the specific framework for fault detection 
and tolerance. The experimental results and discussion are 
presented in Section 4. Finally, the conclusions of the pro- 
posed framework are provided in Section 5. 


Hf. BRIEF DESCRIPTION OF SPND 


An in-core neutron instrument is one of the most important 
types of nuclear power equipment. SPNDs provide a cru- 
cial basis for neutron flux measurement for safe operation, 
treatment of abnormal working conditions, and post-accident 
monitoring of NPPs. In this study, an integrated core instru- 
ment casing assembly was used in a pressurized water reactor 
for third-generation NPPs. As shown in Fig. 1, the core con- 
tains 44 radially distributed measurement channels. In each 
channel, a core detector assembly that is installed in the in- 
strument tube of the fuel assembly contains seven axially dis- 
tributed SPNDs in equidistant layers. 


SPNDs have been developed to meet the requirements of 
small size, long life, and high tolerance to harsh environ- 
ments. The composition of SPNDs is shown in Fig. 2, which 
mainly consists of an emitter, an insulator, and a collector. 
The communication line is constructed using a two-core ca- 
ble in the form of armor and measures the currents generated 
by the emitter and noise currents caused by the background 
core. The currents are then transmitted to the processing cab- 
inets through the connector and transmission cable. Because 
the noise currents are very small, they can be ignored. Gen- 
erally, rhodium SPNDs are used in pressurized water reac- 
tor NPPs. The main generation process of SPND currents 
is determined by the radiation capture reaction, which pro- 
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Fig. 1. Distribution of the neutron flux measurement channels 


ceeds in the emitter material 45 Rh1°3 with the generation and 
subsequent disintegration of the induced beta-active isotope 
45 Rh'°*, Then, the beta particles induced by disintegration 
escape from the emitter with a certain probability and are col- 
lected by the collector, leading to the emitter becoming posi- 
tively charged. Thus, the currents are proportional to the neu- 
tron flux absorbed at the location of the SPND emitter in the 
reactor. 
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Fig. 2. Composition of SPND 


Hl. METHODOLOGY 


The purpose of this study is to use the twin model to en- 
hance SPND detection performance. The main goal of the 
DT model is to create a mirror image of the physical enti- 
ties in the digital realm by observing the data features and in- 
trinsic correlation information among real SPND signals. In 
this study, we conducted a further state analysis of the current 
signals generated by radially distributed SPNDs on the plane 
of the reactor core of a specific NPP. As shown in Fig. 3, n 
SPNDs are uniformly arranged on the cross section of the in- 
core vessel. Each SPND is responsible for the neutron radia- 
tion capture reaction in the relevant space range, which results 
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in n groups of recorded current signals. Each data-twinning 
procedure for an SPND requires signal analysis of a partic- 
ular SPND by observing the signal values of the other n-1 
SPNDs. Therefore, the remaining n-1 SPND signals are used 
as model input characteristic variables, including a particular 
SPND label variable, prior to model training. We can obtain 
n sets of data patterns with related feature and label variables 
for each SPND. Subsequently, a GRNN is employed to train 
and learn the input variables in the data patterns, and the in- 
ternal correlation among the SPNDs can be explained through 
the adaptivity of the neural network. We can obtain n sets of 
sub-models, where the features of the label SPND variable are 
related to the features of the neighboring SPNDs. By integrat- 
ing n sub-models, we obtain an organic whole, which is an 
SPND twin model. Using the twin model, the errors between 
the model outputs and actual signals are calculated to achieve 
SPND fault detection. The model identifies outliers in the 
error sequence in combination with the generalized extreme 
studentized deviate (ESD) statistical test [42], which provides 
rapid SPND fault detection results. A fault tolerance strategy 
should be introduced, the foundation of which is data recov- 


ery for fault variables in response to excessive fault variables 
that result in the failure of fault detection. Fault signals are re- 
placed with normal values through constant troubleshooting, 
data substitution, and data validation, and the performance 
of the twin model is recovered. When all faulty SPND are 
successfully identified, the detection results can assist in di- 
agnostic decisions that can be made to maintain real detector 
assemblies. 


A. Generalized regression neural network 


To deal with sample data comprising complicated multi- 
dimensional variables, prior research concentrated on PCA, 
clustering, and Kalman filter methods to perform signal anal- 
ysis of the SPND unit. These techniques can be used to 
rapidly establish models for axial SPNDs in the same detec- 
tor assembly. In this study, we aimed to analyze the spatial 
correlation between the radially distributed SPND variables 
simultaneously in the plane of the reactor core. The afore- 


mentioned technique is not appropriate for large-scale and 
nonlinear data feature mining; however, the latest neural net- 
work algorithms that an distribute information among neu- 
rons have significant robustness and can quickly analyze com- 
plex nonlinear relationships. As a typical radial basis neural 
network model based on nonlinear regression theory, GRNN 
can process high-dimensional data and mine the mutual ef- 
fects among multiple SPND current features because of its 
capabilities of nonlinear mapping and high learning speed. 
Unlike other types of artificial neural networks, GRNN have 
straightforward structures and training processes. Therefore, 
we used GRNN as the constituent architecture of the SPND 
twin model. As shown in Fig. 4, the GRNN is composed 
of the input, pattern, summation, and output layers. The 
input and output vectors of the corresponding network are 
E Esmee og a and Y = a1, Y2, me n-1 and 
k are the dimensions of X and Y respectively. The number 
of neurons in the input layer is equal to the dimensions of the 
input vector, and each neuron is a simple distribution unit that 
directly transmits the input to the pattern layer. The pattern 
layer is a radial base layer, whose number of neurons is equal 
to the number of learning samples, n-1. The neuron trans- 
fer function P; is typically a Gaussian function and can be 
written as 


P, = exp[—(X — X;)" (X — X;)/207], i = 1,2,...,n— 1 

(1) 
where X; is the training sample corresponding to the 2-th neu- 
ron, T is the transpose of the matrix, and ø is a smoothing 
factor. The smaller ø is, the stronger is the approximation 
ability of P; for the samples. 


There are two types of neuronal calculation formulas in the 
summation layer. One is to sum the outputs of all the neurons 
in the pattern layer. The connection weight between each neu- 
ron and the pattern layer is set to one. The transfer function 
can be expressed as 


S=) P (2) 


The other type of calculation formula is to sum the outputs 
of all neurons in the pattern layer by weight. The transfer 
function can be described as 


Sm; = X wigPi, j=1,2, -k (3) 
t=1 


where w;; represents the weight of the 2-th neuron in the pat- 
tern layer connected to the j-th neuron element in the sum- 
mation layer, S m; represents the summation value of the j-th 
neuron in the summation layer. 


The output layer is calculated by dividing the values ob- 
tained using the two formulas in the summation layer baserd 
on the following equation: 


(4) 


where y; is the output of the j-th output-layer neuron. 
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Fig. 4. Structure of the GRNN 


GRNN can be applied to construct different SPND regres- 
sion submodels that are integrated to form a twin model that 
outputs twin data. 


B. Construction of the twin model 


The state analysis of SPND signals with the application of 
DT can provide a more efficient and intelligent service for 
monitoring the neutron flux distribution and its rate. The 
premise of constructing a twin model is to extract sufficient 
and effective feature information from physical entities. The 
reactor core contains n sets of symmetrically distributed neu- 
tron flux measurement detector assemblies. Although SPNDs 
operate independently of each other, the multidimensional 
currents they produce reveal a strong cooperative link in a 
shared environment. The features of the nearby signals re- 
flect the characteristics of a specific SPND signal. The rela- 
tionship between SPND signals can be explained by the fol- 
lowing equations: 


S= 1 5482, = Sn} 
D = {d1, d2, ...,dn } (5) 
U = {U1 (8), U2(S),.-., Un(S)} 


where S is the set of all SPND individuals at the same height 
in the reactor core, D is the set of sample values for each 
variable in S, and U is the set used to denote the correlations 
between the variables. U; (S) describes how the signal char- 
acteristics of any variable in S are collectively represented 
by the remanent surrounding variables. As seen in Fig. 5, 
the reference network U;(S) is expressed as the connectivity 
between the label variable and other variables. The data of 
the label variable can be twinned by analyzing the network 
of relationships among the variables. By studying the charac- 
teristics of the associated {s1, S2, ..., Sn—1 }, the inherent in- 
formation of the label variables can be obtained. Thus, each 
variable in S can also be collectively described by the oth- 
ers. The twin model is built based on relationship networks, 
and GRNN is used to construct specific submodels to form a 


twin body. By dividing the set D, the n-1-dimensional vari- 
ables are the features J; and the remaining single variable is 
the training label 7;. After an alternate division of n groups, 
the data patterns between J; and T; are expressed as follows. 


[9 də d3 dn—2 dn—1 dn l 
dy d2 d3 : dn— dn dn—1 
dy dy ds ` dn— dn dn—2 
Td =j o. (6) 
dı d3 d4 +++ dn-1 dn | d2 
dz d3 dg --- dn-1 dn | dı 


The division result contains n rows of data patterns, and 
each dimensional variable of SPND signals in set D is des- 
ignated as a label for its corresponding data pattern. For in- 
stance, SPND sn is identified as the label of the first data 
pattern. Similarly, SPND s; is considered as the label of the 
(n+1-i)-th data pattern. To thoroughly explore the correla- 
tion between any target variable and the remaining variables 
in each pattern, GRNN was utilized to train the submodels to 
learn the state relationships between SPND features and la- 
bels. Subsequently, n groups of submodels were built after 
repeated training. An integrated organism M is combined as 
follows: 


M = {M1, Mə, ..., Mn} (7) 


where M; represents the trained GRNN submodels for differ- 
ent data patterns. After learning the algorithm, each GRNN 
submodel outputs the predicted value of the target SPND 
when its input are (n-1)-dimensional feature vectors. All sub- 
models are combined into M, which produces n groups of es- 
timated data and explains the commonality of signal changes 
in all SPND components. The real SPND signals can be 
twinned using these virtual values. The residual between the 
output results of the twin model and the real values can be 
calculated in the given operating scenario to determine the 
faulty SPNDs, and the following status assessment and po- 
tential fault tolerance can then be performed. 
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Fig. 5. The network of relationship U; (S) 


C. Fault detection and tolerance 


Measuring the residual between the outputs of the twin 
model and the actual values allows us to ascertain whether 
faults exist in the detector assemblies. Statistical tests can 
generally be used to determine the deviation points in a resid- 
ual sequence. Traditional statistical techniques are sensitive 
to the presence of outliers because inaccurate data points can 
distort the mean and standard deviation of a data sequence. 
In this section, ESD was employed to detect possible outliers 
in the residual sequence and discriminate the fault from the 
original signals. ESD can maintain a good elimination effect 
when several anomalies exist concurrently in the data. We as- 
sumed that outliers exist in the residual sequence. The maxi- 
mum number of outliers was predetermined as r, and r rounds 
of separate tests were performed by calculating the variable of 
the statistical test, the formula of which is as follows: 


(8) 


where F; is an observed point and E and z denote the sam- 
ple mean and standard deviation, respectively. The observed 
point that maximizes |E; — É| and deviates the most from the 
mean should be located and eliminated during the calculation. 
The aforementioned statistic is then recalculated using the re- 
maining n-1 observed points. The process is repeated until r 
observed points have been removed, resulting in the obtained 
r test statistics R1, Ro, ... , Rr. Corresponding to the r test 
Statistics, r critical values are calculated as follows: 


n—-i+1+t, )(n-i+1) 


pyn—t 


ri = (9) 


( 


where tp» represents the 100p percentage points from the t 
distribution with v degrees of freedom and p = 1 — Snot): 
qa denotes the test confidence level. The number of outliers 
is determined by finding the largest 2 such that R; > A;, and 


thus, the initial hypothesis can be true. 


As mentioned above, the simultaneous failures of exces- 
sive SPNDs in detector assemblies can lead to anomalously 
measured data, which cannot accurately represent the actual 
operational status of the system. A judgment based on fault 
signals will cause the system to fail to reach its normal oper- 
ating condition and even cause personnel to misoperate. The 
capacity of the twin model to perform an effective analysis 
suffers because of the disordered input, and the reliability and 
validity of the detection results cannot be demonstrated by 
the residual information calculated from a broken model. For 
the CNFMS to be dependable and secure, the fault detection 
function must be rehabilitated. In this case, we introduce an 
efficient fault-tolerant strategy, which is implemented before 
multi-fault detection. The most important factor to consider 
when designing a fault tolerance solution is that the infor- 
mation that should be available to the entire system after an 
unexpected failure must be identical to the information that 


would have been obtained if the SPNDs were not faulty. The 
main objective was to achieve data recovery, including trou- 
bleshooting, data substitution, and data verification. 


In the presence of multiple unknown fault sources, fault 
troubleshooting must be performed using multidimensional 
data. A possible selection approach involves selecting the 
vector with the largest calculated residual and performing 
data recovery first. When executing data substitution, the val- 
ues of the faulty variables should be replaced with normal sig- 
nals. To improve the accuracy of the replaced data, the KNN 
classifier was first used to select K instance points from the 
historical database in the nearest neighbor of the feature vari- 
able of the fault points to form a new signal sample. The oper- 
ating condition of K nearest-neighbor points was the closest 
to that of the fault signals, and their variable characteristics 
were in a division range similar to the original signal charac- 
teristics. The weighted average operation was utilized for the 
values of these points, which produced the corrected signal 
values. Finally, the fault signal values were entirely replaced 
by the revised values. 


After data substitution, the Jensen—Shannon (JS) diver- 
gence [43] was employed to confirm the modified data and 
guarantee the consistency of the replacement. It can describe 
the differences between the probability distributions of vari- 
ables, which perform well in the similarity measurement. For 
the two probability distributions, the interval value of the JS 
divergence is [0,1]. The two distributions are more similar 
and vice versa as the value decreases. The degree of JS di- 
vergence can be used to fully evaluate the suitability of the 
modified data. Repeating the preceding steps for a multi-fault 
sample can reduce the number of faulty SPNDs after recov- 
ery until the detection results are satisfactory. When the sat- 
isfied sample is input into the twin model, the output is stable 
and controlled within a small error. In this case, the previous 
fault-detection approach is adequate and fault tolerance is no 
longer required. 


IV. RESULTS AND DISCUSSION 


In this section, the proposed method was evaluated using 
historical monitoring data collected from in-core detector as- 
semblies in a third-generation pressurized water reactor. To 
build a reliable SPND twin model, a sufficient number of 
measurement datasets must be known in advance. The vari- 
ables of the 44-dimensional SPND signals, which were dis- 
tributed radially in a certain layer, were considered as experi- 
mental objects in this study. To improve the generalization of 
the twin model, training data for the twin model were derived 
from the reactor shutdown stage. As shown in Fig. 6, with the 
insertion of control rods, the power of the reactor core grad- 
ually decreased from 100% FP to 0% FP, and all the currents 
of SPNDs decreased. The data for the test set were derived 
from SPNDs during steady operation. 
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Fig. 6. SPND signals during shutdown 


A. Performance of the twin model 


With the continuous reduction in neutron flux during shut- 
down, the number of neutrons captured by the emitter de- 
creased, resulting in a proportional decrease in the currents 
generated by the 44 SPNDs. To extract valuable information 
from physical SPNDs under different power levels and build a 
twin model based on it, submodels for 44 dimensional signal 
variables were executed. First, the training set was segmented 
round after round, and 44 groups of data patterns with several 
feature variables J; and label variable T; were obtained as fol- 
lows: 


dı d2 d3 --- d42 dag | dag 
dy dg dg +++ da2 dag |d43 
dı d2 d3 -> dag daa | dag 
Cy —e « 6a & & lls (10) 
dı d3 d4 --- dag d44| d2 
də d3 d4 --- dag d44| dı 


Each data pattern in (L |T;) was input into GRNN for 
the model to learn the correlation between signal variables, 
and the 44 submodels were well-trained. The superiority of 
GRNN over the other models is ascribed to the convenience of 
its parameter setting. The performance of GRNN can only be 
changed by setting the smoothness factor ø in the kernel func- 
tion. To improve the accuracy of the model, the root mean 
square error (RMSE) was employed as the loss function to 
optimize the value of o for the submodels. Twenty percent of 
the training set was selected as the verification set and a cross- 
verification approach was used to obtain accurate parameter 
optimization results. As shown in Table 1, each submodel ob- 
tained the minimum RMSE with the corresponding setting of 
the optimal o. The optimized parameters helped improve the 
performance of the twin model in the verification set. 


To demonstrate the superiority of GRNN in modeling 
more clearly, we compared the proposed method with sev- 
eral typical approaches, including convolutional neural net- 
works (CNN), multilayer perceptron (MLP), extreme gradi- 
ent boosting (XGB), decision tree (DT), SVM, and Bayesian 
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Fig. 7. Prediction errors of different algorithm-based models 


ridge (BR). All methods were implemented using the same 
data samples from the test set to compare the output accura- 
cies. As shown in Fig. 7, the twin model formed by GRNN 
had the lowest RMSE between the model outputs and actual 
signal values. DT and XGBoost-based models, as common 
tree regression algorithms, had slightly larger errors than the 
GRNN-based model, which can easily overfit SPND signals 
and reduce precision. The learning abilities of the CNN and 
MLP models for SPND signal features were poorer than that 
of the GRNN model. These models are suitable for domi- 
nant feature extraction but tend to ignore the correlation be- 
tween the part and whole. SVR- and BR-based models are 
appropriate for handling small samples instead of analyzing 
larger multi-dimensional signals. Compared to typical ma- 
chine learning algorithms, the twin model established based 
on GRNN can extract rich feature information from isolated 
time-varying data with lower errors and higher prediction 
ability. In the subsequent fault detection phase, the errors be- 
tween the model outputs and actual signals can be used to 
determine the fault location. The estimated values of the twin 
model can help improve the in-core power measurement by 
avoiding interference from fault information. 


B. Fault detection results under different conditions 


The novel contribution of this study is the accurate identifi- 
cation of faulty SPNDs in numerous assemblies using a twin 
model. The signal characteristics analyzed by the twin model 
can help examine the fault characteristics. In this section, 
single- and multi-point SPND fault signals were simulated 
through experiments to determine the effectiveness of fault 
detection under different conditions. The faulty SPND signals 


of diverse measurement channels were simulated by adding 
different percentages of bias to the normal signals based on 
the original operating conditions. 


For single-point fault detection, as shown in Fig. 8, the de- 
viation between the output of the twin model and the actual 
normal value is very small before introducing faults into the 
test set, and the data fit well. In the simulation, a bias of 50% 
of the normal signal value was introduced into SPND#12 at a 
time point after the normal signal. In this case, the twin model 
could maintain a normal output, but the failure of SPND#12 
caused the calculated RMSE to fluctuate significantly at this 
point. The results of the ESD test show that the RMSE of the 
12th point deviates the most, which confirms the existence of 
a faulty SPND#12. 


In reality, complex and changeable interference factors 
usually cause fault signals to vary. Different types of com- 
mon fault behaviors mainly include bias, precision degrada- 
tion, and complete failure, the states of which are abnormal 
owing to the influence of external environmental factors. The 
deviation in the signals also varies at the same power level that 
is used to measure the extent of the magnitude of failures. Ex- 
periments (15 groups x 130 rounds) were conducted to com- 
pare the fault-detection efficiency under different fault types 
and deviation rates. The mean accuracy, an evaluation index, 
is the proportion of the number of successful detections in 
each group of experiments to the total number of experiments. 
Table 2 provides the calculation results for fault detection ac- 
curacy. Except for the detection accuracy of 99.23% in "5% 
drifting", the mean accuracies of other cases were 100%, in- 
dating that the fault variables identified by the ESD test from 
errors are consistent with the actual faulty SPNDs. Therefore, 


Table 1. Parameter optimization results of the 44 submodels 


Submodel o value RMSE (uA) Submodel o value RMSE (pA) 
My 0.3 0.0100 Mi2 0.13 0.0024 
Mo 0.09 0.0020 Mi3 0.07 0.0022 
M3 0.23 0.0069 My4 0.14 0.0029 
M4 0.22 0.0065 Mis 0.33 0.0115 
Ms 0.34 0.0121 Mie 0.02 0.0023 
Me 0.16 0.0038 Miı7 0.42 0.0176 
M7 0.06 0.0023 Mis 0.48 0.0233 
Msg 0.12 0.0021 Mig 0.05 0.0023 
Mo 0.02 0.0023 M20 0.21 0.0061 
Mio 0.44 0.0193 Mai 0.38 0.0146 
Miııi 0.45 0.0202 M22 0.06 0.0023 
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Submodel o value RMSE (uA) Submodel o value RMSE (uA) 
M23 0.32 0.0110 M34 0.05 0.0023 
Moa 0.49 0.0244 M35 0.02 0.0023 
Mo5 0.14 0.0029 M36 0.17 0.0043 
M26 0.44 0.0193 M37 0.07 0.0022 
M27 0.36 0.0133 M38 0.13 0.0024 
Mog 0.35 0.0127 M39 0.05 0.0023 
M29 0.48 0.0233 Mao 0.03 0.0100 
M30 0.25 0.0078 Mat 0.32 0.0110 
M31 0.2 0.0057 Ma2 0.22 0.0065 
M32 0.22 0.0065 M43 0.11 0.0019 
M33 0.31 0.0105 M44 0.14 0.0029 
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Fig. 8. Comparison of the output values of the twin model with the actual values and the RMSE of the models in the two cases 


the proposed method can maintain an extremely high correct 
detection rate for single-point SPND faults under various con- 
ditions. 


Table 2. Fault detection accuracy in various fault conditions 


Mean accuracy with different deviation rates 


Fault types $$ 
5% 20% 40% 60% 80% 
Bias 100 100 100 100 100 
Drifting 99.23 100 100 100 100 
Precision degradation 100 100 100 100 100 


Similarly, it is possible for multiple SPNDs to fail simul- 
taneously. The method for detecting multi-point faults is the 
same as that for detecting single-point faults. However, this is 
still a challenge because multiple faulty signals can contam- 
inate the input of the twin model, leading to rising nondeter- 
minacy of the outputs. In this study, we analyzed the impact 
of the number of faulty SPNDs on the output of a twin model 
and the fault detection accuracy. Because the type and devi- 
ation rate of the failure size had little influence on the detec- 
tion accuracy, the added simulated fault signals were random 


in the multi-point fault experiment (44 groups x 150 rounds). 
The experimental results and comparisons are shown in Fig. 
9. 
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Fig. 9. Influence of the number of faults on the detection accuracy 
of different algorithm-based models 


It was observed that the accuracy of multi-point fault de- 


tection remained at 100% when the number of faulty SPNDs 
was less than 26. In other words, the proposed method can 
perform a relatively efficient detection when approximately 
half of the 44 SPNDs fail. As the number of failures contin- 
ues to increase, the accuracy begins to gradually decrease and 
approaches 50%, which means that approximately half of the 
faulty SPNDs detected by ESD are true. Unacceptable data 
deviating from the normal range is the reason for the occur- 
rence of a sharp fall. Excessive abnormal inputs strongly af- 
fect the original function of the twin model and cause it to col- 
lapse, hindering the realization of subsequent detection. With 
ESD test calculation, the multi-point fault detection accuracy 
of the models based on other algorithms declined earlier than 
that of GRNN-based models under identical conditions. Their 
accuracies decreased and slowed to different degrees with the 
further expansion of failure. Finally, they become more re- 
poseful along their tails. The method proposed in this study 
performs well for both single- and multi-point SPND fault 
detection. 


C. Evaluation of fault tolerance 


The previous analysis proved the superiority of using the 
twin model to detect multiple faults. However, there is still 
a deficiency in dealing with excessive failures, that is, the 
detection accuracy decreases significantly when the number 
of faulty signals exceeds the upper limit. Therefore,efficient 
fault-tolerant solution is required to restore abnormal data. 
For multi-dimensional fault data points, the KNN algorithm 
can help search for K operating points that are closest to the 
faulty signal features in the historical database. On this ba- 
sis, WKNN regression was employed to obtain more accurate 
values by assigning a greater weight to the nearest neighbors. 
The selection of the value of K has a significant impact on the 
results of the regression analysis. This is because an approxi- 
mation error with a smaller K can make the model more com- 
plicated and can likely lead to overfitting, resulting in a larger 
estimation error. According to Fig. 10, WKNN has a lower 
error than KNN without weighting. When K=8, WKNN has 
the smallest RMSE (2.2217uA). Hence, the hyperparameter 
K in WKNN used in this study was set to eight. 


To verify the feasibility of the fault tolerance scheme men- 
tioned above, we selected 100 points of a continuous time se- 
quence under the full-power platform from the test set, and 
random faults were introduced into SPNDs from the 25th 
point to the 100th point in the example. In fault-tolerant pro- 
cessing, data recovery is first performed for the SPND vari- 
able whose RMSE is the largest. Taking SPND#27 as an ex- 
ample, Fig. 11 shows that the deviation between the fault and 
original values is distinct (RMSE=0.0387A), which signifi- 
cantly affects normal system monitoring. From the perspec- 
tive of the twin model, the deviation between the predicted 
and original values was high (RMSE=0.0488WA), making it 
incapable of meeting the requirements of high-precision con- 
trol. In the implementation of WKNN, the deviation distance 
between the regression value and the original values was the 


10 


shortest (RMSE=0.0044,:A), and the data recovery accuracy 
was enhanced compared with that of the twin model. Through 
data recovery, WKNN can perform more reliable correction 
and can maintain stable operation of the system. 


RMSE (pA) 


Fig. 10. Optimization of K in KNN and WKNN 


For data verification, the JS divergence was employed to 
prove the similarity of the probability distribution of the val- 
ues. After the calculation, the divergence values of the fault 
values, predicted values of the twin model, and regression val- 
ues of WKNN with the original values were 4.3409 x 107°, 
7.5058 x10~", and 5.6829x10~", respectively. The regres- 
sion values of WKNN have advantages in terms of similarity 
measurement which verifies the remarkably similar probabil- 
ity distribution and high recovery performance. 
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Fig. 11. Signal value of SPND#27 after data recovery 


We further evaluated the performance of the fault-tolerant 
method on a sample filled with faulty signals. The impact of 
the number of recovered variables on the average RMSE of 
the model outputs and fault detection accuracy is illustrated in 
Fig. 12. In multi-point fault detection, the obtained detection 
accuracy with no measures for recovering data was 52.18%, 
which is the same as before. This is mainly attributed to the 
fact that the output error of the twin model is much larger than 
the normal range, which provides completely false knowl- 
edge and leads to the inability of the ESD test to determine 


that all SPND variables have failed. As the data recovery pro- 
gresses, the faulty variables are substituted by the recovered 
data. Subsequently, the average RMSE decreases, and the 
detection accuracy increases gradually. When the number of 
recovered variables reaches approximately 18, the accuracy 
increases to approximately 100%. In addition, there is no ad- 
ditional need to recover from other failures in this case. 


Detection accuracy (%) 
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Fig. 12. Performance of fault tolerance 


Finally, the effects of the fault tolerance strategy on the 
improvement of twin-model-based fault detection for SPNDs 
are discussed. As shown in Table 3, the number of faulty SP- 
NDs was divided into five intervals in the multi-point fault de- 
tection experiment (150 groups): less than 25, 25-30, 30-35, 
35—40, and more than 40. Without the introduction of fault 
tolerance, the mean fault detection accuracies in the five in- 
tervals were 100%, 72.23%, 53.29%, 52.91% and 53.88%. 
After recovering the faulty data, the number of successfully 
detected faulty SPNDs increased, resulting in accuracies of 
100%. This change confirms that the provision of subse- 
quent fault-tolerant measures to the twin model can efficiently 
enhance the detection accuracy and improve the capacity to 
parse the signal correlation for the model, which is sufficient 
for maintaining healthy monitoring of the system. 


Table 3. Improvement of the fault detection accuracy with fault 
tolerance 
Mean accuracy (%) 


Number of failure 


Before fault tolerance After fault tolerance 


<25 100 100 
25 ~ 30 T223 100 
30 ~ 35 53.29 100 
35 ~ 40 52.91 100 
> 40 53.88 100 
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